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ABSTRACT 

We show that the amplitude of the second acoustic peak in the newly released BOOMERANG-98 and 
MAXIMA-I data is compatible with the standard primordial nucleosynthesis and with the locally broken- 
scale-invariant matter power spectrum suggested by recent measurements of the power spectrum in the 
range 20 — 200 h^ 1 Mpc. If the slope of matter density perturbations on large scales is n » 1, the Hubble 
constant is 0.5 < h < 0.75, and r.m.s. mass fluctuations at 8 h" 1 Mpc are 0.65 < a% < 0.75, then for a 
Universe approximately 14 Gyr old our best fit within the nucleosynthesis bound fish 2 = 0.019 ±0.0024 
requires 0.3 < f2 m < 0.5. Cluster abundances further constraint the matter density to be VL m « 0.3. The 
CMB data alone are not able to determine the detailed form of the matter power spectrum in the range 
0.03 < k < 0.06 h Mpc -1 where deviations from the scale-invariant spectrum are expected to be most 
significant, but they do not contradict the existence of the previously claimed peak at k ~ 0.05 h Mpc -1 , 
and a depression at k ~ 0.035 h Mpc -1 . 

Subject headings: cosmic microwave background - cosmology: theory - cosmology: observations 



1. INTRODUCTION 

The recently released data from BOOMERANG-98 (de 
Bernardis et al. 2000, hereafter BOO) and MAXIMA-I 
(Hanany et al. 2000) have determined the amplitude and 
position of the first and second acoustic (Doppler) peaks 
of the Cosmic Microwave Background (CMB) radiation 
angular spectrum. The first analysis (Balbi et al. 2000; 
Lange et al. 2001; Jaffe et al. 2000) strongly constrained 
the parameter space of cosmological models. The data 
clearly displayed a large first acoustic peak with the max- 
imum at I = l max = 212 ± 7 (Bond et al. 2000) that 
supported the idea of an approximately spatially flat Uni- 
verse with a cosmological constant (vacuum energy) and 
non-baryonic cold dark matter (CDM). However, rather 
unexpectedly, the second acoustic peak at I ~ 500 — 550 
appears to have a low amplitude (especially in the BOO 
data). Another unexpected feature was a shift of the first 
peak to smaller £ than the standard theoretical predic- 
tion lmax — 220 for the flat ACDM model. Even more 
recently, Netterfield et al. (2001, hereafter B01), Lee et 
al. (2001, hereafter M01) have estimated the radiation 
power spectrum up to I ~ 1100 by extending the analy- 
sis of BOOMERANG-98 and MAXIMA-I data to smaller 
angular scales. Halverson et al. (2001, hereafter D01) pre- 
sented the angular power spectrum from the first season 
of DASI observations with data out to I ~ 800. 

A straightforward fit of the earlier data to ACDM mod- 



els with a scale-invariant initial spectrum of adiabatic per- 
turbations lead to: 

a) a high best fit value of the baryon density ilsh 2 » 
0.03 (see, e.g., Tegmark & Zaldarriaga 2000b; Jaffe et 
al. 2000) that significantly exceeds the prediction of the 
standard Big-Bang nucleosynthesis (BBN), updated for 
the recent data on the primordial deuterium abundance: 
n B h 2 « 0.019 ± 0.0024 (Tytler et al. 2000); 

b) a best fit for the density O tot « 1.1 that corresponds 
to a closed Universe (Lange et al. 2001, Jaffe et al. 2000) 
(though the flat Universe is inside the 2a error bars). 

The former result is mainly the consequence of the low 
amplitude of the second acoustic peak, while the latter is 
mostly due to the shift of the first peak to the left. These 
conclusions persist if additional information on the Large 
Scale Structure (LSS) of the Universe, as the density fluc- 
tuation parameter a$, and the shape of the power spec- 
trum of matter are taken into account. In particular, the 
most exhaustive of these efforts made by Tegmark, Zaldar- 
riaga & Hamilton (2000), who performed an 11-parameter 
fit to the current CMB and LSS data, has led to essentially 
the same result about the high baryon density. To avoid 
contradiction with the standard BBN, a number of drastic 
changes in the standard FRW cosmology were proposed, as 
leptonic asymmetry (Lesgourgues & Peloso 2000), delayed 
recombination (Peebles, Seager & Hu 2000), loss of coher- 
ence of primordial perturbations remained from inflation 
(White, Scott & Pierpaoli 2000), admixture of topological 
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defects (Bouchet et al. 2000) or even the absence of dark 
non-baryonic matter (McGaugh 2000). Moreover, even the 
fundamental laws of physics itself, as the constancy of the 
fine structure constant, have been abandoned for the sake 
of explaining these features (Avelino et al. 2000, Battye 
et al. 2001). 

A prcliminar study of the newest D01 data shows no con- 
tradiction with the BBN bound; B01 agrees with the ear- 
lier analysis of Lange et al (2001) except that the baryon 
abundance is reduced to fl^h 2 — 0.027 ± 0.005. This limit 
gets even closer to the BBN estimates if all points are con- 
sidered in the analysis. But at high multipoles, the data 
is not fully consistent. M01 data show a third peak higher 
than the second, while in B01 and D01 both acoustic peaks 
are of similar height. 

All studies quoted above, even those including infor- 
mation on LSS, considered matter power spectra only 
with scale- invariant initial conditions. But Einasto et 
al. (1999c) have shown that the matter power spectrum 
as measured from galaxy and cluster catalogs is inconsis- 
tent with this assumption. A number of different non- 
scale-invariant initial conditions has been recently used to 
analyze the CMB data. First, Kanazawa et al. (2000) con- 
sidered a double inflationary model in supergravity having 
a step in the spectrum located around k = 0.03 h Mpc -1 , 
for which both the matter power at smaller scales and 
the amplitude of the second Doppler peak are reduced. 
A similar matter spectrum was studied in Barriga et al. 
(2000) which resulted from another kind of inflationary 
model with a fast phase transition during inflation. In 
contrast, Griffiths, Silk & Zaroubi (2000) and, most re- 
cently, Hannestad, Hansen & Villante (2000) advocated 
the existence of a bump in the matter spectrum on signif- 
icantly larger scales (k « 0.004 h Mpc -1 ) based on purely 
phenomenological grounds. Einasto (2000) analyzed the 
CMB spectrum resulting from a matter power spectrum 
with a bump at k — 0.05 h Mpc -1 as suggested by Chung 
et al. (2000). 

In most of these papers, the location of non-scale- 
invariant features in the spectrum, or even the functional 
form of these features were introduced ad hoc, with the 
only aim to explain the BOOMERANG-MAXIMA early 
data. In contrast, we adopt a completely different ap- 
proach: without attaching ourselves to any particular the- 
oretical model, we are trying to use previously existing 
observational data as much as possible. To obtain the 
initial (post-inflation) matter spectrum, first, we use the 
empirical present matter power spectrum of Einasto et 
al. (1999a), extracted from galaxy and cluster catalogs 
and estimated in the range ~ 20 — 200 /i -1 Mpc. Sec- 
ondly, on very large scales we assume the initial spectrum 
to be scale- invariant (n = 1) and COBE/DMR normal- 
ized, since this choice produces the best fit to previous 
CMB data for multipoles I < 200 (fortunately, it agrees 
with the prediction of the simplest version of the inflation- 
ary scenario, too). In the next section we shall describe 
two possible ways to match these pieces of the spectrum. 

In this article, we show that such an empirical spec- 
trum gives a possibility to explain the peculiar features 
of the BOOMERANG-98 and MAXIMA-I data without 
changing standard cosmology or the basic laws of physics. 
Also, we do not introduce non-scale invariant features in 
the power spectrum ad hoc, but only those required by 



the observations. In addition to this initial spectrum, we 
assume certain priors: a spatially flat universe, the age 
of the universe between 12 and 14 Gyr and a negligible 
contribution from primordial gravitational waves to the 
COBE/DMR data. We considered models with different 
Hubble constant, cosmological constant, and baryon frac- 
tion. We computed the expected temperature anisotropy 
for each model and found the region of the parameter space 
that best fitted the BOOMERANG-MAXIMA data. In 
Section 2 we describe our methods and results. In Sec- 
tion 3 we find the cluster mass distribution. We calculate 
the initial power spectrum for our models in Section 4 and 
present the main results in Section 5. 

2. TEMPERATURE ANISOTROPIES ON INTERMEDIATE 
ANGULAR SCALES 

Temperature anisotropics are usually expressed in terms 
of spherical harmonics with coefficients ai, n . At a given 
angular scale I, the contribution of the matter power spec- 
trum to the radiation temperature anisotropy is dominated 
by perturbations with the wavenumber k ~ (tt/2)1/Rh, 
where Rh is the radius of the present (post-inflation) par- 
ticle horizon (in particular, Rh = 9900 ft. -1 Mpc for a flat 
Universe with the matter density parameter £l m — 0.3). 
Note that only perturbations with wavenumbers k > 1/ Rh 
contribute to the CMB radiation temperature at any I 
(for I 3> 1 and neglecting the integrated Sachs- Wolfe ef- 
fect which is usually small for scales of interest here). 
Observations in the range I > 200 probe the amplitude 
of the matter power spectrum at recombination on the 
scales also probed by galaxy and cluster catalogs, i.e. on 
scales smaller than 200 hr 1 Mpc. On very large scales 
(k < 0.005 h Mpc -1 ), temperature anisotropies are gen- 
erated by fluctuations of the gravitational potential. The 
amplitude and the slope of the matter power spectrum in 
this region are strongly constrained by the COBE/DMR 
data. To compute CMB temperature anisotropies for 
1 = 2 — 1000, and the matter power spectrum on scales 
10~ 4 < k < 1 h Mpc -1 , we need to provide the initial 
(post-inflation) power spectrum in the whole wavenumber 
range. 

We identify the matter power spectrum with the ob- 
servational spectrum presented in Einasto et al. (1999a) 
(called HD for High Density as the mean power spectrum 
was derived including high-density regions such as clusters 
of galaxies) for k > 0.035 h Mpc -1 . This spectrum is well 
approximated by a power law with the slope n = —1.9 
for wavenumbers larger than k p i = 0.06 h Mpc -1 up to 
k ~ 0.4 h Mpc -1 (the latter value is already sufficient for 
the calculation of CMB multipoles with I < 1000, so we 
use this fit for larger values of k, too): 

P HD (k) = P pt (k/k^y 1 - 9 h- 3 Mpc 3 , (1) 

where P p i = 1.02 x 10 4 (cr 8 /0.65) 2 h -3 Mpc 3 is the value 
of the power spectrum at k = k p i (see Table 2, Einasto 
et al. 1999a). In the range 0.035 < k < 0.3 h Mpc -1 
this spectrum differs from the observed power spectrum 
of galaxies, by a scale factor or bias. On shorter scales, 
k > 0.3 h Mpc -1 , the observed spectrum is corrected 
for non-linear evolution effects (for details see Einasto et 
al. 1999c); these corrections do not introduce any sig- 
nificant differences to the CMB spectra at the scales of 
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interest. Note that in this range, the spectrum of Einasto 
et al. (1999a, 1999c) is in good agreement with results 
obtained by other authors, for recent new data see Miller 
et al. (2001), among others. Since we do not know the 
exact value of the bias between the galaxy and matter 
power spectra, the value of k p i is a free parameter of our 
model which is expressed in terms of as once the shape 
of P(k) is fixed. On the other hand, on very large scales, 
k < k m — 0.03 h Mpc" 1 , we accept the scale-invariant 
(n — 1) COBE normalized spectrum since it gives the 
best fit to CMB data for I < 200. The location of k m is 
chosen to be in agreement with the observational cluster 
spectrum of Miller & Batuski (2000) and the de-correlated 
IRAS Point Source redshift catalog (PSCz) galaxy spec- 
trum of Hamilton, Tegmark & Padmanabhan (2000). 

In the intermediate range, k m < k < k p i, error bars 
are large and there is no complete agreement between dif- 
ferent authors about the exact form of P(k). For this 
reason, we use two different models of P(k) in this re- 
gion: one more conservative and another based on the 
Einasto et al. (1999a) data. In the first case we use at 
small scales the power law behavior given in ([l]). This 
spectrum is extrapolated until it crosses the linear COBE- 
normalized spectrum that is based on the primordial scale 
invariant spectral index n = 1 and evolved through recom- 
bination using cosmological parameters discussed above. 
On large scales the spectrum based on is n = 1 used. 
It was calculated using the CMBFAST program. This 
power spectrum is called HDL (L standing for linear ex- 
trapolation). In the second case we apply the observa- 
tional HD spectrum from Einasto et al. (1999a) up to 
k = 0.035 h Mpc -1 , i.e., up to the smallest wavenumber 
for which it contains reliable data. Then P(k) is linearly 
matched to the scale invariant n = 1 COBE normalized 
spectrum at k = k m = 0.03 h Mpc -1 . For smaller k, we 
use the scale invariant spectrum as in the first case. The 
resulting spectrum has a distinct bump at k = kb, and 
a depression in the whole region k > k m , as compared to 
the COBE normalized n = 1 spectrum, see Fig la,b. Thus 
we call it HDB (B for bump). Notice that the difference 
between the two spectra is still within la error bars. 

The HDL spectrum is similar to those considered in 
Kanazawa et al. (2000) and Barriga et al. (2000). The 
HDB spectrum is more complicated, it is not monotonous. 
Its form is in agreement with the power spectrum pro- 
posed by Einasto et al. (1997); the existence of the 
bump at k = 0.05 h Mpc -1 (now understood to be a 
relative bump above a depression in the spectrum) could 
help to explain some quasi-periodic features in the large- 
scale distribution of Abell clusters noted in the latter and 
other papers. It is interesting that the double feature - 
the bump at k = 0.05 h Mpc" 1 and the depression at 
k ~ 0.035 h Mpc -1 - are seen in the Miller & Batuski 

(2000) and Hamilton, Tegmark & Padmanabhan (2000) 
data, too, though within la error bars. Similar features 
have been reported in two recent preprints: Hoyle et al. 

(2001) found a small bump at k ~ 0.06 h Mpc" 1 (assuming 
fi TO = 0.3) in the power spectrum of the 2dF QSO Redshift 
Survey, and Silberman et al. (2001) claimed a wiggle in the 
power spectrum with an excess at k ~ 0.05 h Mpc" 1 and 
a deficiency at k ~ 0.1 h Mpc" 1 , based on peculiar veloci- 
ties of galaxies from the Mark III and SFI catalogs. Bumps 
at k w 0.05 h Mpc" 1 and k ~ 0.2 h Mpc -1 and depres- 



sions (valleys) at k ~ 0.035 h Mpc" 1 and k ~ 0.1 h Mpc" 1 
have been confirmed by Miller, Nichol & Batuski (2001) in 
their analysis of power spectra of Abell and APM clusters 
of galaxies, and PSCz galaxies. Note, however, that there 
are no such features in the recent REFLEX X-ray cluster 
survey (Schuecker et al. 2000). Certainly, both HDL and 
HDB are non-scale-invariant spectra. Though we have ar- 
rived at them using purely empirical arguments, we shall 
discuss how they could have arisen in inflationary models 
in Section 4. 

To compute the radiation and matter power spectrum 
at the present epoch we used the CMBFAST program de- 
veloped by Seljak & Zaldarriaga (1996). Given a set of pa- 
rameters and scale-invariant initial conditions of the type 
5(k) = Ah^-V/ 2 for the rms Fourier amplitudes of the 
density contrast at horizon crossing, we obtained a matter 
power spectrum with spectral index n at large scales. Pa- 
rameter A was normalized to reproduce the COBE/DMR 
amplitude. In this work, we assumed that a difference 
between the power spectrum predicted by a given cosmo- 
logical model today, P(k), and the observed power spec- 
trum, Pffo(fc), was due to initial conditions only. Then, 
as all scales of interest are still in the linear regime, the ini- 
tial rms Fourier amplitude of matter density perturbations 
that give rise to the observed spectrum is 

S HD (k) = d(k)(P HD (k)/P(k)y/ 2 . (2) 

In this expression, the power spectra are evaluated at the 
present time, and the amplitude of SjjD(k) is evaluated at 
the horizon crossing. 



Our main goal is to determine if the B01 and M01 
new measurements agree with the observed matter power 
spectrum and with the standard BBN. Therefore, we re- 
stricted our analysis of the parameter space to a rather 
limited set, centered upon the region that best fitted 
the CMB data before BOOMERANG-98 and MAXIMA-I 
data (see Tegmark & Zaldarriaga 2000a). We considered 
spatially flat models with and without cosmological con- 
stant (fi m + Qa = 1, Q m = Qcdai + ^b), with matter 
density in the range 0.3 < Q m < 1, with n — 1 and baryon 
fraction 0.005 < 57s/i 2 < 0.035 centered on the range sug- 
gested by Tytler et al. (2000). Particular attention was 
paid to models with a cosmological constant. We took the 
age of the Universe to to be between 12 and 14 Gyr, in 
agreement with the recalibration of the cosmic ages made 
by Feast & Catchpole (1997) using new Hipparcos distance 
determinations, and from estimations of the age of globu- 
lar clusters (Jimenez 1999). The resulting Hubble constant 
varied from h ~ 0.5 to 0.75, consistent with observations. 

Models with a small mixture of massive neutrinos have 
also been discussed in the literature, but this component 
has little effect on the radiation power spectrum in the 
range considered here. However, it is important to realize 
that as neutrino free-streaming erases the matter power 
spectrum at small scales, they have an effect on the ini- 
tial matter power spectrum computed using (||). If the 
dark matter is "cold" and Q m ~ 1 then the initial spec- 
trum deviates strongly from scale invariant initial condi- 
tions. This is expected since the standard CDM has too 
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Fig. 1. — (a,b) Matter power spectrum of ACDM models with scale invariant (dashed line), HDB (solid line) and HDL 
(thick dot-dashed line) initial conditions. The matter density is Q m — 0.5 in (a) and 0.3 in (b). The box encloses the 
region where the observed power spectrum is extended to match smoothly the scale invariant spectrum at large scales. 
Diamonds and bars represent the observed matter power spectrum Pnoik) and the associated la error bar (data taken 
from Table 2 of Einasto et al 1999a). This area has been enlarged in (c,d) to show the difference between the spectra 
HDL and HDB. For clarity the scale invariant ACDM model is not plotted. (e,f) Radiation power spectrum of the same 
models. Filled diamonds and open circles correspond to B01 and M01 data, respectively. 

Table 1 

Position of the first maximum of the CMB spectrum 
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Q m = 0.3 


Scale invariant 
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217 


HDL 
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HDB 
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much power on small scales. If we assume that the ob- served spectrum deviates only weakly from the Harrison- 
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Fig. 2. — Goodness-of-fit contours of the likelihood function at 95% (thick solid line) and 99% (thin solid lines) confidence 
levels. The fit was performed using the newly released data B01, M01 up to I ~ 650. No calibration uncertainties were 
included. Contours on the left side of panels correspond to scale invariant ACDM models and are shown for comparison; 
right contours correspond to HDL or HDB models. All models correspond to a Universe age of 14 Gyr and to matter 
power spectra with a n = 1 spectral index on large scales. Upper two panels: solid lines correspond to a HDL power 
spectrum with (a) as = 0.65 and (b) <j% = 0.75. Lower two panels: solid lines correspond to a HDB power spectrum with 
(c) er 8 = 0.65 and (d) as = 0.75. Dot-dashed straight lines lines correspond to the la BBN limits on the baryon fraction. 



Zeldovich spectrum, then the theoretical power spectrum 
on scales smaller than the bump (k > 0.05 h Mpc -1 ) co- 
incides with the observed spectrum only for a restricted 
set of cosmological parameters. In particular, demanding 
that the observed spectrum deviates only slightly from 
scale invariance, we fix the matter density parameter to 
be fi m w 0.3 (if the non-baryonic dark matter is "cold"). 
A small mixture of "hot" dark matter would allow for a 
larger matter density (Einasto 2000 and Section 4 below). 

In Figure 1, we present the COBE-normalized matter 
power spectra (upper panels a and b with an enlarged ver- 
sion in the middle panels c and d) and radiation power 
spectra (lower panels e and f) for three different types of 
spectra: scale invariant (dashed line), HDB (solid) and 
HDL (thick dot-dashed), for a universe with to = 14 Gyr, 
the baryon fraction Sigh 2 = 0.02 and the matter den- 
sity f2 m = 0.5 (left panels) and f2 m = 0.3 (right pan- 
els). In this figure we use as — 0.75 in the left panel and 
as = 0.65 in the right panel. The HDL power spectrum 
follows Eq (0) till k ~ 0.08 h Mpc -1 (for VL m = 0.5), and 
till k ~ 0.06 (for fi m = 0.3) where it follows the ACDM 
spectrum. Then HDL coincides or is below the scale in- 



variant power spectrum. The HDB spectrum follows the 
observed spectrum up to k — 0.035 h Mpc -1 , and then 
is linearly extrapolated until it reaches the scale invariant 
spectrum at k — 0.03 h Mpc -1 . It can be above or below 
the scale invariant spectrum. In the radiation spectra, the 
main differences between HDB and HDL are around the 
first Doppler peak. The data on the matter power spec- 
trum was taken from Einasto et al. (1999a), while the 
CMB data is the latest M01 and B01 (this is also the data 
used to fit the models shown in Figure 2). 

From the initial power spectrum, we compute the radi- 
ation power spectrum and compare the expected level of 
anisotropy with the CMB data. We used the radical com- 
pression of the cosmic microwave background data method 
described in Bond, Jaffe & Knox (2000). We computed 
the parameter space compatible with the newly released 
B01, M01 and, for comparison, with the earlier data recal- 
ibrated according to Jaffe et al. (2000). We restricted our 
study to I < 650 since it allows a better comparison with 
earlier parameter estimates. Further, above that scale the 
error bars are most affected by ±13% effective beam uncer- 
tainty and this uncertainty is not included in the published 
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Fig. 3. — Cluster mass functions for models with the density parameter fi m = 0.3 and fl m — 0.5 are plotted in the left and 
right panels, respectively. The observed cluster mass functions are given according to Bahcall & Cen (1993) (thick solid 
line indicates the data and thin solid lines the width of the error corridor) and Girardi et al. (1998). Model functions are 
given for scale-invariant CDM models (with a cosmological constant), and for our models HDB with a bump for a 8 = 0.65 
and CTg = 0.75. 



error bars. In Fig 2 we give the contours at the 95% (2a) 
and 99% confidence levels for ACDM models with scale 
invariant, HDB and HDL power spectra. The figure cor- 
responds to a fit to the B01, M01 data. On the upper axis 
we plot the corresponding Hubble constant. All models 
presented in Fig. 2 correspond to the slope n = 1 on very 
large scales. Figure 2a shows that scale invariant models 
are marginally consistent with the BBN bounds only at the 
99% level. Our results also indicate that HDL and HDB 
models with to = 12 Gyr (not shown in the figure) require 
£1a > 0.75 (at the 2a level), hardly compatible with direct 
measurements of Sl m = 1 — f^A- On the contrary, models 
with to = 14 Gyr lay inside the standard BBN limits and 
fit the recent CMB data very well for both the HDL and 
HDB spectra. We repeated the same analysis with the 
original recalibrated data. In this case, the contours were 
wider and scale invariant models become compatible with 
the BBN bound. Let us remark that, as seen in Figure 1, 
BOOMERANG-98 and MAXIMA-I have different ampli- 
tudes at the scale of the second acoustic peak, thus it is 
harder for the models to fit them both. However, in both 
cases the contours for HDL and HDB spectra were cen- 
tered on the same region of the parameter space, implying 
that our results are robust. 

In Table 1 we give the position of the maximum of the 
first acoustic peak for the models presented in Fig. 1. Let 
us remark that since, in general, the HDB and HDL spec- 
tra have less power than a scale invariant model with the 
same cosmological parameters, the position of the maxi- 
mum of the first acoustic peak is shifted with respect to the 
scale invariant model. This effect is easily explained by the 
k — I correspondence given at the beginning of this section: 
the depression in the HDB spectrum relative to the scale- 
invariant spectrum in the range k = 0.03 — 0.05 h Mpc -1 
results in the decrease of the corresponding angular radia- 



tion spectrum for Z = 200 — 500 (if Vl m = 0.3). This effect 
explains the position of the first Doppler peak without in- 
voking a small positive spatial curvature. 

3. CLUSTER MASS FUNCTION 

Another important constraint of cosmological models is 
the cluster mass function. This function is rather sensi- 
tive to cosmological parameters and features in the power 
spectrum on intermediate and small scales which have the 
highest weight in the cluster formation process. To check 
our results we calculated the cluster mass function using 
the Press-Schechter (1974) method. Mass functions were 
found for two sets of models; we used the Hubble con- 
stant h = 0.65 and the baryon fraction Jig = 0.05, the 
matter density was taken Q m = 0.3 and fl m — 0.5, vary- 
ing the density of the cold dark matter and the vacuum 
energy (cosmological constant) respectively. The Hubble 
constant used is a compromise between recent new esti- 
mates by different teams (Parodi et al. 2000, Sakai et 
al. 2000). As before, we used only spatially flat models 
(f2 m + Q\ = 1); the slope of the spectrum on small scales 
was taken to be n = —1.9, and the spectrum amplitude 
parameters as — 0.65, 0.75 were used. 

The results of calculations are shown in Figure 3. The 
HDB model with Q m — 0.3 and as — 0.65, that fits the 
CMB data rather well, has a too low abundance of clusters 
while the scale invariant ACDM model and HDB model 
with Q m = 0.3 and as — 0.75 lie within the range allowed 
by observations. It is not surprising that the last two mod- 
els coincide since the ACDM model with Q m = 0.3 has 
£78 = 0.78, very close to our HDB model with Q, m = 0.3 
and as = 0.75. This as value is in good agreement with 
other independent estimates (Einasto et al. 1999b, Van 
Waerbeke et al. 2001). The ACDM model with O ro = 0.5 
has too high abundance of clusters over the whole mass 
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range (about a factor of ten). In contrast, Q m = 0.5 HDB 
models with both values of <7g have cluster mass functions 
which lie within the allowed region. In other words, cluster 
mass functions of HDB models are rather insensitive to the 
density parameter. This result is expected since our HDB 
models with both density values (but an identical Hub- 
ble parameter) are identical on medium and small scales, 
k > 0.06, which have the highest weight in determining 
the cluster mass function. 

Using the preprint version of this paper, Gramann & 
Hiitsi (2001) investigated properties of models similar to 
our models HDL, fixing the break of the power spectrum 
at k p i = 0.05 h Mpc -1 (instead of at 0.06 h Mpc -1 ). 
They confirmed our results that the amplitude of the sec- 
ond and third peak in the radiation power spectrum de- 
creases and fits the CMB data better. Gramann & Hiitsi 
also calculated the cluster mass function with the Press- 
Schechter method and found that the amplitude of the 
power spectrum of their models was too low which lead 
to a low cluster mass function incompatible with observa- 
tions. They suggested to use an initial power spectrum 
which has a depression around k = 0.1 h Mpc -1 and an 
increased amplitude starting from k = 0.2 h Mpc -1 . Mod- 
els with such spectra were investigated by Suhhonenko 
& Gramann (1999). Such spectra also yield satisfactory 
agreement both with the observed cluster mass function 
and with the recent CMB spectrum data. A similar power 
spectrum was found in the recent paper by Silberman et 
al. (2001) based on the analysis of galaxy peculiar veloci- 
ties, and by Miller, Nichol & Batuski (2001) based on the 
analysis of power spectra of Abell and APM clusters of 
galaxies, and IRAS Point Source redshift catalog galaxies. 

To bring our HDB model with fl m = 0.3 and erg — 0.75 
into a better agreement with the CMB data, perhaps a 
model with a second feature on small scales is needed, as 
suggested by Gramann & Hiitsi (2001) and Miller, Nichol 
& Batuski (2001). The latter authors considered the pres- 
ence of two bumps and two depressions (valleys) in the 
observed power spectrum as evidence for baryonic fluc- 
tuations, as discussed by Eisenstein et al. (1998). How- 
ever, this hypothesis encounters one problem: the posi- 
tion of the first baryonic bump on the matter power spec- 
trum is expected at the wavenumber k m 0.07 h Mpc -1 , 
whereas the observed bump is located at the wavenum- 
ber k « 0.05 h Mpc -1 . The bump at k = 0.05 h Mpc -1 
corresponds to the scale of the supercluster-void network, 
~ 120 /i -1 Mpc, as seen in the distribution of high-density 
regions in a pencil-beam near galactic poles (Broadhurst 
et al. 1990), and in the distribution of clusters of galax- 
ies (Einasto et al. 1997a, 1997b). More accurate data 
are needed to decide whether the observed features of 
the power spectrum are due to a high baryon fraction 
or are primordial effects. At present, we are inclined to 
accept the second alternative since, without requiring new 
physics, it fits the CMB data within the BBN bounds quite 
naturally. 

4. POST-INFLATION INITIAL CONDITIONS 

Though we did not require a priori that our initial spec- 
tra should be "explicable" by any theoretical model, it ap- 
pears that, for Q m rj 0.3, the form of the HDB spectrum 
suggests an inflationary model with a fast phase transi- 
tion occurring approximately 50 e-folds before the end of 



inflation, similar to those with a sudden jump in an infla- 
ton potential which were considered in Starobinsky (1992) 
and Adams, Ross & Sarkar (1997); see also the recent 
paper by Adams, Cresswell & Easther (2001) (which ap- 
peared after the preprint version of the present paper was 
sent to the archive) where a very similar spectrum is ob- 
tained. As for the simpler HDL spectrum, and for the 
HDB spectrum with Q m w 0.5, it is rather typical for dou- 
ble inflationary models, see, e.g., the models in Polarski 
& Starobinsky (1992), Gottlober, Miicket & Starobinsky 
(1994), as well as the models proposed in Kanazawa et al. 
(2000) and Barriga et al. (2000). Also, a similar step- like 
spectrum appears in the case when the inflaton potential 
has a jump of its first derivative (Starobinsky 1992). Thus, 
relatively small changes in cosmological parameters and in 
the present form of P(k) may lead to significantly different 
models of the phase transition during inflation. 

To elaborate this point further, let us compute the ratio 
of power spectra: 

S{k) = P HDB (k)/P CDM (k) . (3) 

Here, we use the power spectrum with a bump, Phdb {k) , 
and the scale- invariant ACDM power spectrum PcDAi(k), 
calculated for the same set of cosmological parameters as 
the spectrum with a bump. The function S(k) charac- 
terizes the deviation of our accepted spectrum from the 
scale-invariant case, namely, the initial spectrum Po(k) oc 
kS(k). In Fig. 4 we plot the ratio S(k) for a Hubble con- 
stant h = 0.65, baryon density FIb = 0.05, and erg = 0.65 
and erg = 0.75 for three different values of the matter den- 
sity. We do not plot the ratio using the HDL spectrum 
since S(k) is the same as for HDB for k > 0.06 h Mpc -1 
and coincides with the scale invariant case (S(k) — 1) for 

k ~ 0.05 h Mpc -1 as explained in Section 2. 

By construction, on large scales the HDB spectrum 
is identical to the scale- invariant ACDM spectrum, i.e., 
S(k) — 1. Deviations start at k ~ 0.03 h Mpc -1 : all initial 
spectra have a depression around k — 0.035 h Mpc -1 , a 
bump at k = 0.05 h Mpc -1 , and an approximate power law 
dependence on smaller scales. The initial spectrum given 
in (||) depends only slightly on the er 8 parameter. For the 
matter density f2 m = 0.3, deviations of our accepted spec- 
tra from the scale-invariant spectrum are moderate and lo- 
calized near k = 0.05 h Mpc -1 . We call these spectra "lo- 
cally broken-scale-invariant". In contrast, for the matter 
density fl m = 0.5 and fl m = 1.0, deviations on small scales 
become large. Even if we do not include a small fraction of 
hot dark matter, the function S(k) can be brought closer 
to unity on the scale interval 0.06 < k < 0.4 h Mpc -1 
using a tilted spectrum. Tilted spectra have been often 
used to explain the observed power spectrum of galaxies. 
This helps to bring S(k) rj 1 over the whole range only 
for fi m = 0.3. For higher values of the matter density, the 
tilt needed becomes too large and is excluded using the 
COBE data, that suggest n = 1 ± 0.1. 

To summarize, the CMB data in combination with the 
matter spectrum and cluster abundance favor the density 
parameter Q m ~ 0.3. For this density value, the initial 
matter power spectrum is approximately scale invariant. 
Figure 2 indicates that for Q, m ~ 0.3 the CMB data prefers 
(78 = 0.65 for both the HDL and HDB spectra, while the 
value of fi m ~ 0.5 is preferred by the data when erg = 0.75. 
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h=0.65, n b =0.05 h=0.65, a b =0.05 
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Fig. 4.— The initial power spectrum, expressed as the ratio of our adopted HDB spectrum to the scale-invariant ACDM 
spectrum, S(k) = PHDB(k)/PcDM{k). The initial power spectrum is calculated for 3 values of the matter density, 
fl m = 0.3, 0.5, 1.0. In the left panel erg = 0.65, in the right panel er 8 = 0.75. 



However, this density value exceeds considerably the value 
preferred by other data sets (Perlmutter et al. 1998, Riess 
et al. 1998, see also Ostriker & Steinhardt 1995), and also 
requires large deviations from the scale-invariant initial 
conditions. 

5. CONCLUSIONS 

We have shown that the amplitude of the second acous- 
tic peak measured by BOOMERANG-98 and MAXIMA-I 
is compatible both with the standard BBN and with the 
matter power spectrum obtained from galaxy and cluster 
catalogs, if we accept a non-scale-invariant power spec- 
trum. If we assume that the age of the Universe is 14 Gyr, 
and the initial power spectrum has an index close to the 
Harrison- Zeldovich value, n = 1 ± 0.1, then our best fit 
within the standard BBN bound gives for the matter den- 
sity parameter, 0.3 < VL m < 0.5. If we demand that the 
matter power spectrum yields a correct cluster abundance, 
then the matter density is further constraint: fi m « 0.3. 

The amplitude of the second peak in the radiation power 
spectrum and the shape of the matter power spectrum can 
be both simultaneously explained using models with the 
standard baryon abundance. Since our accepted spectra 
have less power than the scale-invariant model with the 
same cosmological parameters (and before the new data 
was available) we had predicted a low third acoustic peak. 
This testable prediction of our models can be clearly seen 
in Figs 1(e), 1(f): the third acoustic peak is not higher 
(and typically lower) than the second peak. While the 
original paper was being refereed, Netterfield et al. (2001) 
extended the analysis of the BOOMERANG-98 data to 
higher multipoles and indeed found an amplitude of the 
third Doppler peak being of a similar amplitude as the 
second one. However, Lee et al. (2001) claim to have seen 
a larger third peak in the MAXIMA-I data. In this re- 
spect, let us remark that HDL and HDB spectra coincide 
for k > 0.05 h Mpc -1 , their CMB angular spectra are prac- 



tically the same for I > 500 (if Q TO = 0.3) while in scale 
invariant models with a high baryon abundance and mod- 
els with leptonic asymmetry the third peak is higher than 
the second. As the data around the third peak improves, 
its amplitude will be an important test on our models. 

Our second result is the shift of the first acoustic peak 
to a smaller value of I as compared to the scale-invariant 
spectrum (see Table 1). This effect is especially notice- 
able in the case of the HDB spectrum. This brings the 
BOOMERANG-MAXIMA data to excellent agreement 
with a spatially flat Universe, removing the need for a 
slightly positive spatial curvature. Also, the cluster mass 
function for the HDB type models is well within the ob- 
servational error bars. 

Finally, it is important to notice that the data do not 
contradict the existence of a bump and a depression of 
the HDB spectrum. The existing CMB data are not suffi- 
cient to determine the detailed form of the matter power 
spectrum in the range 0.03 < k < 0.06ft. Mpc -1 where we 
expect a deviation from the scale-invariant spectrum to 
be most significant. Analysis of the data around the third 
peak will be useful to determine the amplitude of the mat- 
ter power spectrum at small scales, but is not relevat for 
the question about the feature. But before this analysis 
can be carried out, the discrepancy between Maxima and 
other experiments should be resolved in some way, i.e., is 
the third peak high or low. 
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